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We study the excitation of spatial patterns by resonant, multi-frequency forcing in systems un- 
dergoing a Hopf bifurcation to spatially homogeneous oscillations. Using weakly nonlinear analysis 
we show that for small amplitudes only stripe or hexagon patterns are linearly stable, whereas 
square patterns and patterns involving more than three modes are unstable. In the case of hexagon 
patterns up- and down-hexagons can be simultaneously stable. The third-order, weakly nonlinear 
analysis predicts stable square patterns and super-hexagons for larger amplitudes. Direct simula- 
tions show, however, that in this regime the third-order weakly nonlinear analysis is insufficient, 
and these patterns are, in fact unstable. 
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The variety of patterns and planforms that have been 
observed in surface waves on vertically vibrated fluid sur- 
faces (Faraday waves) is remarkable [11, [3, d, . As elu- 
cidated in various theoretical investigations [J [1,01 this 
extreme variability is a result of the fact that the wave 
form of the vibration's forcing function allows for a de- 
tailed tuning of various aspects of the interaction between 
different plane-wave modes, which can stabilize complex 
patterns like super-lattice patterns and quasi-patterns. 
Motivated by this richness of patterns we are investigat- 
ing here the effect of time-periodic forcing with different 
wave forms on systems undergoing a Hopf bifurcation to 
spatially homogeneous oscillations. 

In order to describe a forced Hopf bifurcation within 
a weakly nonlinear framework the forcing must be suffi- 
ciently weak. For it to have a qualitative effect on the 
system it must then include frequencies that are close 
to one or more of the low-order resonances of the sys- 
tem, i.e. its spectrum has to contain frequencies close to 
the Hopf frequency ujh itself (l:l-forcing), close to 2L0h 
(2:l-forcing), or close to 'iwh (3:l-forcing). These strong 
resonances lead to additional terms in the weakly non- 
linear description and qualitatively affect the system [1]. 
To wit, in the weakly nonlinear regime the complex oscil- 
lation amplitude A satisfies a complex Ginzburg-Landau 
equation of the form 



dA 

— = [pi + iv)A^(\ + ili)l\A-(\ + ia)A\A\^ + 



'-fA + r]A^ + C- 



(1) 



The forcing terms C,, jA, and rjA^ express the effect of 
forcing the system at the frequencies w, 2w, and Sw, re- 
spectively, with = w/i + I . The parameter /x expresses 
the distance from the Hopf bifurcation, which is shifted 
by a Od??!^) compared to the unforced case. Here we 
will focus on the case C = 0. To include the forcing near 
the 1 : 1-resonance one can eliminate the inhomogcneous 
term C by shifting A by the fixed-point solution Aq sat- 
isfying C -(m + 11^)^0 + (1 + ia)Ao\Ao\^ - 7A0 - rjAl 
and use Aq instead of C as an external parameter [^ . 



As is apparent from ([T]) , the forced Hopf bifurcation is 
described by an equation that is very similar to a two- 
component reaction-diffusion equation. The only and 
significant difference is the term involving f3 which char- 
acterizes the dispersion of unforced traveling wave so- 
lutions, which would be absent in the reaction-diffusion 
context. It plays, however, an essential role in excit- 
ing patterns with a characteristic wavenumber [l^ and 
cannot be omitted (cf. ^ below). Pattern selection in 
a general two-component reaction-diffusion system has 
been studied in detail by Judd and Silber [ll[, who find 
that in principle not only stripe and hexagon patterns 
can be stable in such systems, but also super-square and 
super-hexagon patterns. They also find that despite the 
large number of parameters characterizing these systems 
surprisingly only few, very special combinations of the 
parameters enter the equations determining the pattern 
selection. 

Amplitude Equations. In this paper we will stay 
below the Hopf bifurcation taking ^ < 0. Thus, as in 
Faraday systems, in the absence of forcing, no oscillations 
arise. To investigate the weakly nonlinear stable stand- 
ing wave patterns possible in ([T]) we derive amplitude 
equations for spatially periodic planforms. The linear 
stability of the state A = is easily obtained by splitting 
the equation and the amplitude A into real and imagi- 
nary parts {A = Ar + iAi). The usual Fourier ansatz 
Ar^i oc e*'^^ yields then the neutral stability curve jn{k) 
with the basic state being unstable for 7 > 7„(fc). The 
minimum 7c (fc) of the neutral curve is found to be at 



,2 _ fi + i^P 



7c 



(1+/32) 



(2) 



Since /j, < 0, the condition > implies that spatial 
patterns arise only if the detuning of the forcing relative 
to the Hopf frequency is such that waves with non-zero 
k are closer to resonance than homogeneous oscillations 
with fc = fld\. A typical neutral curve is illustrated in 
Figure[T]for ^ = -l,i^ = 4, /3 = 3 and C = 0. The weakly 
nonlinear analysis presented in this paper is valid for val- 



2 



ues of 7 near 7c. The range of validity is restricted by 
7„(fc = 0) where spatiaUy homogeneous oscillations are 
excited by the forcing, which interact with the standing- 
wave modes with wavenumber kr. 




Figure 1: Neutral stability curve for ([1} with = — 1, /3 = 3, 

i/ = 4 and C = 0. The critical point (fcc,7c) ~ (y*^' 'Tu)) '® 
marked by a red circle. 

To determine the stability of the various planforms we 
first determine the amplitude equations for rectangle pat- 
terns, which are comprised of two modes separated by an 
angle 9 in Fourier space. We expand (A^, A^) as 

(3) 

where < e <C 1 and the complex amplitudes Zi (T) and 
Zg{T) depend on the slow time T = e^t. The wavevectors 
are given by ki = (fcc,0) and kg — {kcCoa{9),kcSm{6)). 
We also expand 7 as 7 = 7c + £^72 . 

The usual expansion leads to the amplitude equations 
for {Zi,Ze), 



dZi 
'dT 

dZg 

It 



= x{-/-jc)Zi-{bo\Zi\^ + h{e)\Zg\^)Zi, (4) 

- Xi-f^Jc)Zg-{bii9)\Zi\^ + bo\Zg\'')Zg. (5) 



If = ri G Z, the quadratic nonlinearity induces 
a secular term and the expansion has to include three 
modes rotate by 120° relative to each other. 



A, 



V2 



+ 0(e2). 



(6) 

The parameters can be chosen such that the quadratic 
solvability condition is delayed to cubic order, yielding 



dZi 
'dT 



Ml ~ lc)Zi + (jZ2Z:i - 
-(6o|^iP + 62(1^21' + 1^31'))^! 
with similar equations for Z2 and Z3. 



(7) 



More complex patterns can be described by combining 
these two analyses. For example, a super-hexagon pat- 
tern comprised of two hexagon patterns {Zi, Z2, Z^} and 
{Z^, Z^, Zq} that are rotated relative to each other by an 
angle 9sh is described by the amplitude equation 



dZi 



=A(7 - jc)Zi + aZ^Z-j,- 
-(6o|^iP + 62(|^2p + |^3|'))^i 



J=0 



and corresponding equations for Zj, j = 2, . . . , 6. 

The coefficients of the amplitude equations can be writ- 
ten in a simple form, setting rj = rjr + irif. 



A 



271+^(077, + 77,) 



b,{9) = 6V + 8/(^^)x, 
62 = 6^Ij + 1Qx + 4 



with 



-2Q(Q-/3)v/r+^ 

/3(l + a2) 



(8) 

(9) 
(10) 
(11) 

(12) 
(13) 
(14) 



^ -Pjiy - M/3) 2 
^ 2(M + ^^/?)''' 

4(1 + (3^)^, [{a + /?) r?, - (a/3 - 1) rj, 
/32(l + a2)(z.-^/3) 

and the angle dependence given by 

^ 3 + 16cos4g 
' ~ (4cos2 0-l)2- 

Here a = \/\ + /32 +/3. The scaling of the nonlinear coef- 
ficients is based on a normalized eigenvector v — 

Pattern Selection. As shown by Judd and Silber [11| 
for general two-component reaction-diffusion systems, at 
the point of degeneracy at which the quadratic coefficient 
a and with it the coefficient x vanishes not only stripe 
patterns but also hexagon or triangle patterns can be 
stable. The conditions for hexagons (or triangles) to be 
stable are 




< 

20 ^ , ^ 

< w < . 

15 ^ 3 



(15) 
(16) 



Whether hexagon or triangle patterns are stable depends 
on higher-order terms in the amplitude equations p^ . 
which are not considered here. For < V < ~T5' 
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hexagons are unstable to stripes. Both patterns bifur- 
cate unstably for < 0. 

Comparing the conditions (|15I16|) with expressions 
(|12|13|14p shows that in the system considered here stripe 
patterns and hexagons can be stable at the point of de- 
generacy, depending on the system parameters a and f3. 
Specifically, at the degeneracy one has rji = —arjr, imply- 



ing 



(a = 0) = - 



4a??^(l + /3^)^ 



< 0. 



(17) 



Here we have made use of the fact that (3{i' — fijS) > 0, 
which follows from < and the condition > 0. 
Whether condition is then satisfied depends on the 
specific unforced system: only a certain range of values of 
a leads to stable hexagons, which is determined by the 
a— dependence of V'- A distinguishing feature of these 
hexagon patterns is that both 'up'- and 'down'-hexagons 
are simultaneously stable and are likely to form compet- 
ing domains. Figl^] shows an example of the competition 
between 'up'- and 'down'-hexagons in a numerical simu- 
lation of (P). 

Unfolding the degeneracy, i.e. taking < \a\ 1, 
the transition to hexagons becomes transcritical and 
hexagons are stable to stripes for a 7— range given by 



4A(6o + 262 



< 7 < 7c + 



a^(2bo + 62 )_ 
A(62 - boy 



iHS- 



Note that ^hs > 7c, even if (|15ll6p are not satisfied, 
since stripes do not exist for 7 < 7c. The instability 
of hexagons at 7^5 only arises if &2 > ^Oi that is, if 
3ip+lAx/^+4> > 0. With (T 7^ the up-down symmetry of 
the amplitude equations ([7]) is broken and, depending on 
the sign of a, either up- or down-hexagons are preferred. 

Turning to other planforms, Judd and Silber found 
that rectangular planforms cannot be stable at or near 
the degeneracy point [ll|- Interestingly, however, they 
find that while super-hexagons cannot be stable at the 
degeneracy point, they can arise in a very small param- 
eter regime in its vicinity if the conditions 



> 0, 

— < tp < — 
21 ^ 3 



(18) 
(19) 



are met. They then can be bistable with hexagons. We 
find that in our system the condition ^TE\i cannot be met. 
Rewriting (j) in terms of the deviation a from the degen- 
eracy condition. 



aP{v-fiP) V 



(g -I- I3)r]i 



shows that - for small |(t| - can be made positive only 
by making rji small as well {rji = 0{a)). Even then cj) 
can only be slightly positive, (f> = 0(cr), requiring that 
-0 = 0{a) in order to satisfy Under these conditions 



all cubic coefficients in (I9I10I11|) would become of 0{a) 
and without knowledge of the next-order coefhcients no 
stability predictions can be made. 

Away from the degeneracy, for a = 0(1), the above ar- 
guments suggest that it should be possible to satisfy the 
stability conditions (|18ll9p . For a = 0(1) they are, how- 
ever, not the correct stability conditions since they ignore 
the angle dependence of the cubic coefhcients, which is 
0{a). We use (I18I19P therefore only as a guide to locate 
parameter regimes in which super-hexagons may be ex- 
pected to be stable and then determine the full stability 
eigenvalues. One such case is given by the linear pa- 
rameters used in Figure [T] with the nonlinear parameters 
a = —1 and rj ~ e'^/"*. 

Moreover, for a = 0(1) the weakly nonlinear analysis 
predicts also stable rectangle patterns satisfying |6i(6')| < 
bo if 



a > 



1 + / 76 , , 



Within the hexagon sub-space, for 1 62/^^0 1 > 1 hexagons 
become unstable to stripes as 7 is increased beyond jhs- 
Satisfied simultaneously, the two conditions I62/&0I > 
1 and |foi(0)| < bo would yield a situation in which 
hexagons are unstable to stripes, which in turn arc un- 
stable to rectangle patterns. 

Using direct numerical simulations of the forced com- 
plex Ginzburg-Landau equation ([1]) we have studied to 
what extent the predictions of the weakly nonlinear anal- 
ysis are borne out. In the degenerate case cr = we find, 
as predicted, either stripes or hexagons to be stable de- 
pending on the value of a chosen. A typical hexagonal 
pattern obtained from random initial condition is shown 
in Fig|21 As expected, it exhibits competing domains of 
up- and down-hexagons. 




Figure 2: Competing domains of up- and down-hexagon do- 
mains starting from random initial conditions with linear pa- 
rameters as in Fig. [T] and nonlinear parameters a = —1, 
rir — 0.4 and rji — 0.4(v'T0 + 3). rjr and r)i are chosen so that 
a = 0. 
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Figure 3: Comparing amplitude from numerical simulation 
with amplitude predicted by weakly nonlinear analysis. Sim- 
ulation parameters: = = 4, /3 = 3, a = —1, rj = e^^^*. 
FPn corresponds to the fixed point obtained from numerical 
simulation, FPa to that obtained from the weakly nonlinear 
calculation to cubic order. 

Away from the degeneracy, a = 0(1), the validity of 
the weakly nonlinear analysis can be severely restricted 
by the fact that the amplitudes of all stable branches are 
0(1), which formally suggests the significance of higher- 
order terms in the expansion. Indeed, in and near the 
parameter regimes for which the weakly nonlinear analy- 
sis predicts stable super-hexagon patterns we do not find 
any indication of their stability. To assess explicitly the 
significance of the higher-order terms in the amplitude 
equations for a = 0(1) we extract them directly from 
the numerical simulation for the case of hexagon pat- 
terns. Fig. [3] shows the numerically determined depen- 
dence of \A\~'^ d\A\/dt on \A\ for 7 = 7c. For very small 
1^1 it agrees well with the result obtained by the third- 
order weakly nonlinear theory, which yields the straight 
dashed line. However, even for 7 = 7c the amplitude \A\ 
saturates only at a value of \ A\ « 0.14 for which the third- 
order theory deviates significantly from the full result. A 
fit of I (i| A|/c?t to a higher-order polynomial shows 
that the magnitude of the quintic term in the amplitude 
equation reaches a value of 20% of the cubic term. The 
fixed point obtained from weakly nonlinear analysis to 
cubic truncation deviates from the numercially obtained 
fixed point by 30%. This supports out interpretation 
that in this regime the cubic amplitude equation does 



not allow quantitative predictions. 

In summary, we have investigated the regular spatial 
planforms that can be stably excited in a sytem under- 
going a Hopf bifurcation by applying a periodic forcing 
function that resonates with the second and third har- 
monic of the Hopf frequency. We have done so within 
the weakly nonlinear regime by deriving the appropri- 
ate amplitude equations describing the selection between 
various planforms. By tuning the phase of the forcing 
close to 3w/i one can always reach the point of degen- 
eracy where no quadratic terms arise in the amplitude 
equations, despite the quadratic interaction in the un- 
derlying complex Ginzburg-Landau equation. Hexagons 
arise then in a supercritical pitchfork bifurcation. De- 
pending on the parameters of the unforced system we 
find that in this regime either the hexagon patterns or 
stripe patterns can be stable. In the former case com- 
peting domains of up- and down-hexagons are found in 
numerical simulations when starting from random initial 
conditions. Surprisingly, despite the extensive control 
afforded by the two forcing terms, no square, rectangle, 
or super-hexagon patterns are stable in the vicinity of 
this degeneracy, irrespective of the parameters of the un- 
forced system. Only in the regime in which hexagons 
arise in a strongly transcritical bifurcation the weakly 
nonlinear theory predicts the possibility of stable rectan- 
gles or super-hexagons. There, however, direct numer- 
ical simulations of the complex Ginzburg-Landau equa- 
tion indicate no such stability and we show that terms of 
higher order in the amplitudes are relevant. 

By introducing a further forcing frequency, which is 
also close to the 2:l-resonance, the transcritical bifur- 
cation to hexagons can be avoided. As we show in a 
separate publication, the corresponding, more elaborate 
weakly nonlinear theory then correctly predicts stable 
quasi-patterns comprised of four, five, or more modes 

For larger values of the forcing the spatially periodic 
standing wave modes interact with a spatially homoge- 
neous oscillation that is also excited by the forcing. The 
interaction between these two types of modes could lead 
to interesting patterns, which are, however, beyond the 
scope of this paper. 
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